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Abstract. Supernovae remnant shock waves could be at the origin of cosmic rays up to energies in excess of the 
knee (_E ~ 3 ■ 10^^ eV) if the magnetic field is efficiently ampli fied by the streaming of accelerated particles in the 
shock precursor. This paper follows up on a previous paper jPelletier et al lliood) which derived the properties 
of the MHD turbulence so generated, in particular its anisotropic character, its amplitude and its spectrum. In 
the present paper, we calculate the diffusion coefficients, also accounting for compression through the shock, and 
show that the predicted three-dimensional turbulence spectrum k±S3d{k\\, k±) oc k'^^kj" (with fcy and k± the 
wavenumber components along and perpendicular to the shock normal) generally leads to Bohm diffusion in the 
parallel direction. However, if the anisotropy is constrained by a relation of the form fcy oc k^'^ , which arises when 
the turbulent energy cascade occurs at a constant rate independent of scale, then the diffusion coefficient loses its 
Bohm scaling and scales as in isotropic Kolmogorov turbulence. We show that these diffusion coefficients allow to 
account for X-ray observations of supernova remnants. This paper also calculates the modification of the Fermi 
cycle due to the energy lost by cosmic rays in generating upstream turbulence and the concomittant steepening 
of the energy spectrum. Finally we confirm that cosmic rays can produced an amplified turbulence in young SNr 
during their free expansion phase such that the maximal energy is close to the knee and the spectral index is close 
to 2.3 in the warm phase of the interstellar medium. 

Key words. Physical processes:Acceleration of particles-shock waves-turbulence-Interstellar medium: Supernova 
remnants 



1. Introduction 

Supernovae (SN) blastwaves are probable sites of cosmic 
ray (CR) ac celeration up to energies of a few lO^"* eV as 
discussed bv lLaeage fc Cesarskv I (|1983j) . The cosmic rays 
are accelerated through the diffusive shock acceleration 
mechanism which invokes upstream and downstream tur- 



bulence and whose efficiency is fully determined by the 
transport prope rties of t he co s mic rays in these region s 
[see reviews by IPrury I lll983|l. Ijones fc El ison I lll99l|) , 
and more recently bv lDrurv I l)200l[l . iHillas 1 1)2005|) ]. The 

maximum cosmic-ray energy attainable in supernova rem- 
nants (SNR) scales as the product of magnetic field and 
size of the accelerator; the above Lagage & Cesarsky es- 
timate for the maximum energy assumes that the flue- 
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tuating part of the magnetic field 5B has saturated at a 
level of the order of the fa r upstream interstellar value Bao 
i McKenzie fc Volk Ill982t) . lowever, recent observational 
and theoretical studies have indicated that young SNR are 
likely sites of efficient magnetic field amplification with a 
fluctuating component reaching levels ^ -Boo- 
High resolution X-ray observations of young SNR have 
revealed the presence of sharp external rims associated 
with the blastwave, pr oduced by synchrotron radiation 
of relativistic electrons llpotthelf et al.ll200lHHwang et al.1 



l2002t iLong et al.ll2nO,'^Tvinkll20n4t lBalletll200,'^^ As dis 

cussed in Section 2 below, the smallness of the rim size im- 
plies a downstream diffusion coeffici ent for the relativistic 
electrons close to its B o hm value |Berezhko et al.lE^003l 
iBerezhko fc V olk" "2004; ' Parizot et al I '2006') i.e £>„ ~ 
rr c/3. lBell fc Lucek. (,200lfl and lLucek fc Bell . (,2001.) have 
investigated in some details the non-linear development 
of the cosmic ray streaming instability and obtained an 
amplified magnetic field well in excess of Boo- These au- 
thors have suggested that the tangled character of the 
magnetic field implies a CR mean free path of the order 
of its Larmor radius (the Bohm diffusion regime). Both 
Bohm diffusion and a strong magnetic field amplification 
seem necessary to push the maximum cosmic-ray energies 
up to the CR knee at a few PeV or even up to t he CR 
ankle at a few EeV l|Ptuskin fc Zirakshvilh lf2n0,'i 120051) . 
This amplification mechanism then appears essential to 
our understanding of galactic cosmic-ray phenomenology 
and theory. 

However, previous theoretical work suffers from var- 
ious limitations. For instance, the streaming of cosmic- 
rays amplifies forward waves at the same time as it damps 
backward waves, and the transfer of energy between these 
two ty pes of waves has not been modeled bv lBell fc Lu^ 
l|20ni[) . The numerical simulations using coupled "Particle 
in Cells" and magnetohydrodynamics codes suffer from 
limited wave number dynamics and a crude modeling of 
the acceleration mechanism, questioning the value ob- 
tained for the non-linear saturation level of the magnetic 
field. The argument in favor of a Bohm diffusion regime 
is also phenomenological since no MH D theory has rigor- 
ously predicted such regime yet fs ee ICasse et al. ( 200^ 
and references therein], although iLucek fc Bell I |200l() 
did observe isotropisation of the CR distribution on a 
timescale of a gyroperiod in their numerical experiments. 
Fi nally, the non-resonant instability mechanism uncovered 
has not been taken into account in previ- 
ous studies. These points have motivated us to investigate 
in some details the generation of turbulence in the shock 
precursor by cosmic-ray streaming through both resonant 
and non-resonant instabilities (PcUcticr et al 2006), here- 
after refered as paper I. We have notably found that the 
three-dimensional spectrum of upstream turbulence so ex- 
cited is anisotropic and scales as Szd{k\\,k±) oc fc|~^A:J"~^ 
(with fc|| along and perpendicular to the shock normal 
and to the mean magnetic field) , and that the downstream 
turbulence is further amplified by the jump through the 



shock. 

In the present work, we propose to use the new results of 
paper I to investigate the acceleration of cosmic rays in 
young SNR. In particular, we explicitly calculate the dif- 
fusion coefficients both upstream and downstream using 
Monte-Carlo simulations and analytical quasi-linear the- 
ory. We thus provide a connection between the prediction 
of turbulence spectra and the observations of electron syn- 
chrotron radiation. We also apply the results to the Fermi 
cycle of cosmic rays. Finally we follow in detail the energy 
budget, accounting in particular for the energy spent in 
amplifying the upstream magnetic field. 

2. Observations 

2.1. X-ray and radio observations of ttie outer rims 

Quite recently the X-ray instrument Chandra has imaged 
the continuum emission (4 to 6 keV) in different yo ung su- 
pernova remnants li ke Cas A ^Cotthelf et a,lJ2nmf) . Tycho 
ijHwang et al]l2002t) . The X-ray radiation originates from 
very thin sheets behind the blast wave. The same con- 
clusion have been drawn by XMM-Newton ob servations 
of the Kepler SNr l)Cassam-ChenaY et al. lf2004|) . The de- 
rived physical width of the emission region is less than 
4" or 2 10^^ cm in Cas A, for a shock velocity V^h — 
5000 km/s, less than 4" or 1.5 10^^ cm in Tycho for 
a shock velocity P^h — 4600 km/s. The results for three 
young SNr (Cas A, Tycho and Kepler) and two old er ones 
(SN1006 and G347.3-0.5) have been presented in iBallefi 
|2005.) . The effective rim size ma y however be overest i- 
mated due t o projection effects IjBerezhko et al.1 
iBalletl l|2005f) has also shown that this radiation is in- 
consistent with thermal models in a uniform medium 
which predict emission everywhere up to the interface, 
with only a slight maximum at the blast wave. A second 
argument against thermal models is the rather feature- 
less spectrum in the 1-10 keV band which appears to be 
inconsistent with the X-ray brightness [see for instance 
ICassam-Chenai et al. I ()2004f) in the case of the Kepler 
SNr] . The only other possible source of X-ray radiation is 
suprathermal particles. Non-thermal bremsstrahlung (by 
low-energy suprathermal electrons) could still be compat- 
ible with the observed rims for a low thermal gas temper- 
atur e of ^ 1 ~ 2 10^ K and a high enough target den- 
sity l)Balletll2005j) . Another natural mechanism is the syn- 
chrotron radiation by u ltra-relativistic electrons with ener - 
gies of few tens of TeV ((Vink fc Lam ing'200a':'Vink"2004l). 
The radio synchrotron emission from GeV electrons is ex- 
pected to have larger extent as the loss timescale at these 
energies is much longer. T his seems to be ind eed the case 
for C as A llCotthelf et a,1.lEo01.l and SN1006 ijl.ong et a,lJ 
120031) (for this source the rim profile has been plotted in 
radio and in two X-ray bands showing a systematic sharp- 
ening when shifting to the highest frequencies). The case 
of Tycho SNr seems to be more complex, here the two rims 
(radio a nd X-rays) do n ot track each other in brightness 
well (iHwang et al.l We shall discuss in details the 
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physical consequences of these observations on electrons 
and cosmic rays in section 

2.1.1. Constraints from the shock precursor 

The shock precursor region is of prime importance as it is 
the place where mos t of the turbulence is expected to be 
generated. lAchtcrbc rg et alJ (^994) noticed that an upper 
limit on the precursor size implies an upper limit on the 
electron mean free path and a lower limit on the level of 
turbulence. Depending on the scale over which the emis- 
sivity drops upstream ^1/2, the electron mean free path 
is Ae < 3 {Vsh_/c) i!i/2/[acos^(0_Bn)] where Osn is the 
obliquity angle between the shock normal and the up- 
stream magnetic field, and a accounts for an emissivity 
drop; it equals l/log(2) in an exponential profile or 3 for 
a smoother profile where the turbulence is self-generated. 

With the help of hi gh resolution radio observations, 
lAchterberg et all l)l994|) concluded to a magnetic field 
amplification of about one order of magnitude com- 
pared to the standard interstellar values. This conclu- 
sion is however not correct in case of non diffusive trans- 
port where supplementary information about th e num- 
ber of upstream shock crossings are required l|R,agotl 
If the magnetic field is purely turbulent upstream, 
one can replace in the previous expression cos^6'b„ by 
1/3. This leads to upper limits for Ae of the ordrer 
of 10-3pc (V;h/3000 km/s) (^i/a/O.Olpc). 



A. 



For instance, 
< 1.1 10-3- 



pc, 

pc, Ae < IQ-^pc 



Tycho and Kepler have £1/2 — 7 10 pc 
and l^h ^ 5400 km/s, l-^/^ ~ 7 10"^ 
respectively. The diffusion coefficient for radio electrons 
D ~ 'i;Ae/3 is limited to 3 1025-26 cm^/s. The Bohm 
value L»Bohm ^ 6 10^^ (E/lGeV) (B/l/iG)"^ cmVsevenif 
compatible with these constraints only provides loose con- 
straints. In X-rays the emitted radiation must be asym- 
metrical due to the magnetic field jump at the shock 
front. The observed X-ray profile is strongly dependent 
on the projection effects and it appears that most (if not 
all) of the X-ray phot ons should come from downstream 
l|Berezhko et al.ll2nn3j) . 

To conclude, it is important to note that the synchrotron 
interpretation of the radio and X-ray emission is a nat- 
ural option to explain the high resolution images from 
X-rays and radio telescopes and the non-thermal spec- 
tra by Chandra and XMM-Newton. But this model still 
have some difficulties and should be considered with cau- 
tion. However, we adopt in the following the synchrotron 
mechanism as the dominant radiative process in the sharp 
outer rims. 



2.2. Consequences on the particle acceleration process 
2.2.1. Electrons 

The size of the X-ray rims downstream the blast wave have 
important constraints on the relativistic particle trans- 
port. For electrons, one can expect that in the loss limited 
case the size of the rim is set fixed either by advection 



or by diffusion. In fact, the electron energy is expected to 
be close to the maximum energy given by tcooi(i'max) — 
iacc(-E'max) — '^/^sh' thdA, casc the advcctiou scale 
^adv = icool Vsh should bc close to the diffusive scale 
£riiff = \/k tcooi l)Berezhko et al.l l2003l IVink fc Lamind 
\i E ~ i?max, by comparing the advection length 
with the rim sizes inferred from the observations, one de- 
duces typical magnetic field strengths of the ord er of a 
few 1 00 /iG in most of all young SNR cited above l|Balleli 
l200fil) . pointing towards an important magnetic field am- 
plification at the shock precursor by the accelerated par- 
ticles themselves. Assuming a diffusion coefficient D = 
a Dsohm =i a X 6 10^4 (E/10 TeV) (B/lOO /iG)-i cmVs, 
with a > 1, the previous co nstraints impose a between 1 
and 10 |Parizot et alll2006(l . 



2.2.2. Cosmic-rays 



In the previous paragraph, we have found that high energy 
electrons accelerated close to the Bohm value produce syn- 
chrotron X-ray profiles compatible with the most recent 
Chandra observations of young SNr. Likewise, protons and 
heavier nuclei are prone to diffusive shock acceleration in 
this highly disordered magnetic field. However, one may 
expect protons energies to be much higher than few tens 
of TeV and have Larmor radii allowing them to explore 
the largest turbulent scales. At these wavelengths no infor- 
mation is available on the characteristics of the turbulence 
yet, i.e. the spectral index, the turbulence level. Even if the 
diffusion coefficient is close to its Bohm value at the en- 
ergy where the electrons produce synchrotron photons in 
the 4-6 keV range, there is no observational evidence that 
it remains tr ue at higher energies . In most of the previous 
approaches jBell k. Lucekl I2OOI : IPtuskin &: Zirakshvilli I 



l2003t iBerezhko et al.ll2003t 16611120041) only an heuristic 
argument was used, i.e. the particle mean free path is ex- 
pected to match the Larmor radius in a completely disor- 
dered magnetic field. Numerical experiments have however 
shown that in the limit of a strongly turbulent magnetic 
field, in the case of a Kolmogorov turbulence, the diffusion 
coefficient approaches its Bohm value at large rigidities 
only, i.e. fo r Larmor radius cl o se to the max i mum turbu- 
ill l2002j) . IParizot et all ||200(tI) using 



lence scale 

the X-ray observations from a sample of young SNR have 
investigated all possibilities testing the spatial transport 
of cosmic-rays in an isotropic turbulence. They conclude 
that spectral turbulence indexes smaller than 3/2 (index 
of the Id spectrum) are to be rejected as they lead to 
a diffusion coefficient at the maximum proton energies 
smaller than the Bohm diffusion coefficient. Even if the 
diffusion regime is close to Bohm, the authors did not 
find maximum cosmic-ray energies beyond the CR knee 
(~ 3 lO-'^^eV), questioning the capability of young SNr to 
produce the high energy part of the galactic CR spectrum. 
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3. Summary of the results of paper I 

In order to discriminate between the different transport 
regimes described above, we have investigated in paper 
I in greater detail the turbulence spectrum that results 
from cosmic-ray interactions in the upstream medium. 
This work has shown that the turbulence growth oc- 
curs as a result of two instabilities. At large dist ances 
from the shock front, the non-resonant instability |BeLl 
l2004|) that results from the non-zero return current in the 
thermal plasma, amplifies the turbulence at short wave- 
lengths (as compared to the typical cosmic-ray Larmor 
radius at that distance from the shock front). At closer 
distances, cosmic-rays of Larmor radius tl excite res- 
onantly forward turbulent modes of wavenumber fcy = 
l/r-L ©ell k Lucck 200 IJ. These instabihties govern the 
growth of the parallel forward modes and the concomit- 
tant damping of the backward parallel spectrum. While 
the non-resonant growth saturates when energy is redis- 
tributed through non-linear transfers as fast it is input, 
the resonant instability is quenched by advection. The 
non-resonant instability that leads to turbulent spectrum 
5(^11 ) oc acts at larger distances from the shock 

front than the resonant instability. This latter then su- 
persedes the former, leading to a turbulent spectrum at 
the shock front 5(^11 ) oc This one-dimensional spec- 

trum is related to the three-dimensional spectrum by 
5(^11) = (27r)~^ J d^fc_LS'3d(fc||, /c^), and is normalized ac- 
cording to / dA:||S'(fc||) — 2ttSB'^ /B"^, where SB and Boc 
denote respectively the turbulent component and the ISM 
magnetic field. 

On a timescale that is shorter than the instability 
growth time and the advection time, energy is distributed 
by resonant three-wave interactions in the perpendicu- 
lar direction k±. This redistribution leads to anisotropic 
three-dimensional turbulent spectra fcj.S'sd oc fcJ"/C||~^, 
with /3 = 1 as before, and a = {7 - 2(3)/?, = 5/3 if 
the energy transf er rate of the cascade is a constant inde- 
pen dent of s cale llGoldreich fc Sridharll 19951: iGaltier et al.1 
I2OOO; Galtier et al. Il2005|) . This also implies a relation- 
ship between perpendicular and longitudinal wavenum- 
bers: k^^L^^-^ oc (fcj^imax)^^'^, where Lmax is the maximum 
turbulence scale. These results have bee n found in agree- 
ment with recent numerical simu lations ijCho fc Vishniaci 
l2nO(lHMaron fc Cxoldreich ''200lV 

It has also been shown in paper I that interactions between 
two Alfvcn and one slow magnetosonic waves provide ef- 
ficient transfer of energy along the longitudinal direction, 
which does not arise when only interactions between three 
Alfven waves are considered. This allows to maintain the 
backward parallel spectrum close to the level of the ampli- 
fied forward spectrum. It has also been pointed out that 
the non-resonant instability induces a left-right symmetry 
breaking which provide suitable ground for further ampli- 
fication by a dynamo. This aspect however has not been 
studied in detail and is left for further work. 
The turbulence that is generated by cosmic ray inter- 
actions is compressed along the shock normal at shock 



crossing. In particular, the perpendicular magnetic field 
components are amplified by the shock compression ratio 
r while the longitudinal component remains unchanged. 
Correspondingly, the turbulent modes wavenumbers 
are increased by the same factor r, which corresponds to a 
similar decrease of the turbulence coherence length in that 
direction, while perpendicular modes are unchanged. The 
downstream turbulence is thus even more anisotropic. In 
the following, we discuss the implications of these turbu- 
lent spectra with respect to cosmic ray transport at the 
shock front, hence with respect to Fermi acceleration. 

4. Diffusion coefficients 

The diffusion coefficients due to resonant pitch angle scat- 
tering of cosmic rays with turbulent modes can only be 
derived analytically in the limit of weak turbulence (quasi- 
linear theory, Jokipii 1966). Since one expects a high tur- 
bulence level in the vicinity of the shock front, the diffusion 
coefficients have to be computed irumerically using M onte- 
Carlo techniques as described in lCasse et al.l ((20021) . The 
main results of this paper are that the pitch angle fre- 
quency follows the same scaling as the quasi-linear theory 
(hereafter QLT) namely i^s/wl ~ rjp^~^ for p < 1, where 
(3 is the index of the turbulence spectrum, — 1 /ts is the 
scattering frequency (and Tg the scattering time), and cjl 
the Larmor frequency. For p > 1 the ratio I's/'-^L decays 
as 1/p. The reduced rigidity is defined as: 

P = 't^^^ = fcininT-L , (1) 

with i,jiax the maximum length of the turbulence and /cmin 
the associated minimum wavenumber. Thus, the parallel 
diffusion coefficient _D|| ~ (1/3)w^Ts oc e^'^ for small p 
(e is the particle energy) and oc for large p, v c 
is the cosmic-ray velocity. The transverse diffusion coeffi- 
cient D± has an unusual scaling which neither corresponds 
to that from quasilinear theory or to the Bohm coefficient, 
as it is rather controlled by field line chaos, leading to 
= ?7^^^D|| (with e ^ 0.3). In numerical experiments, 
there is no observed Bohm scaling, only a Bohm maximum 
for the scattering frequency when p ^ 0, 1 and 77 ^ 1. 

These scaling laws are modified when anisotropic tur- 
bulence is generated, as will be seen in the next section. 
We have performed new numerical simulations of parti- 
cle diffusion in turbulent magnetic fields, paying particu- 
lar attention to the anisotropic nature of the turbulence 
expected in the vicinity of the shock front. We have de- 
composed the three-dimensional turbulence on separate 
discrete grids of perpendicular and parallel wavenumbers. 
This and the required number of Monte-Carlo runs result 
in rather long computation times. We have also derived 
the diffusion coefficients using the quasi-linear theory per- 
turbative methods and used these to interpret the numer- 
ical diffusion coefficients, as discussed below. 

Assuming the existence of a mean field component {B) , 
one defines the pitch angle diffusion frequency Vs — = 
{Aa'^)/At, where a is the pitch angle with respect to the 
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mean field {B) and Aa the random jump of a during At, 
with Tc <C Ai <C Ts . The correlation time Tc is of order the 
Larmor time and the latter strong inequalities hold only 
in weak turbulence regime or in strong turbulence regime 
for p ^ 1. The function g = VsI^'l can be expressed in 
terms of the thre e-dimensional turbulent spectrum [see 
ICasse et alJ ()2002l) and appendix A] as: 

^S3^m,k^)R{k,p) , (2) 

where Ax(t) and Lu{k) stand for the variation of the un- 
perturbed particle trajectory within a timescale r and for 
the mode pulsation respectively. R{k,p) is the resonance 
function defined by: 



R{k,p) 



dKr) e[»'^-'^"-^'"('=)-"l cosKr) 



(3) 



A detailed derivation of the function g including estimates 
in the stro ng turbulenc e regime can be found in the ap- 
pendix of ICasse et alJ (|2002) . The Landau-synchrotron 
resonances in Eq.© are obtained by inserting the gyro- 
motion, and we get 



R{k,p) 



E 

n— — oo 



J^{k±ri^ sin a) x 



{S [fc||rL/i - (71 + l)]+S [fc||rLM +{n- 1)] }4) 

The expression obtained for linearly polarized Alfven 
waves in Eq.dll ) is identical to oth er derivations 
l)Schlickeiseill200i lYan fc Lazarian Il2002|) . 



4.1. Upstream d if Fusion coefficients 
4.1.1. Modified weak turbulence spectrum 

We start with a model of weak MHD turbulence 
l|Galtier et alJl2000() upstream modified by the streaming 
instability of cosmic-rays. In that case (see paper I), the 
3D turbulence spectrum reads: 



S3dik\\,k±_) = 5*0 Un ^ 



where M|| = k\\ £\\ and u± = kj 



(5) 

. Out of generality we 
leave the parallel and perpendicular i± magnetic field 
maximal lengths unspecified; to make contact with pre- 
vious notations, in isotropic turbulence they would corre- 
spond to the maximal length Lmax introduced earlier. For 
a given Larmor radius tl we have thus two rigidities p|| 
and p±. 

Since we are mostly interested in the scaling of the 
diffusion coefhcients with rigidity or Larmor radius, we 
drop the numerical prefactors for convenience. Inserting 
Eq.p into Eq.gl) we find 



g cx V |n± 1| 



Only the harmonic n — contributes to the integral sig- 
nificantly. The function g is then found to be constant for 



poo 

L 









• • • • • 
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Fig. 1. Parallel (solid dots) and perpendicular (diamonds) 
diffusion coefRcients in units of tlc as a function of rigidity 
27rrL/imax, Lmax denoting the largest scale of the turbu- 
lence. The turbulence level is 77 = 0.95 and the turbulent 

5/3 1 

spectrum k±S3d oc fc|| . The Bohm diffusion regime 
is clearly apparent for the parallel diffusion coefficient. The 
dotted lines have slope ^1/3, and 1. On can see a slight 
offset in the ordinates of the data points between the two 
regimes p < 10~^ and p > 10~^; this offset is a numerical 
artifact, see text for details. 



tl <C i±, i.e. the pitch-angle scattering time is propor- 
tional to the particle energy. In other words, a turbulence 
spectrum S oc l/fc|| as generated by cosmic-ray interac- 
tions in the shock precursor leads to Bohm diffusion with 

This result agrees quite well with the numerical com- 
putation of the diffusion coefficient in the high turbu- 
lence level limit, shown in Fig. ^ for rj — 0.95, where 
T] = SB^/{SB^ + {B)"^). As the turbulence level decreases, 
the scaling remains the same but the numeri cal prefactor 
for the diffusion coefRcients increases, see also lCasse et alJ 
(2002). This figure reveals a slight offset in the ordinates 
of the data points between the two regimes p < 10^^ and 
p > 10~^, which is a numerical artifact. In order to ob- 
tain data points down to p 10~^ when only a limited 
dynamic range in wavenumber is allowed by computing 
limitations, we have conducted two separate simulations, 
one for each rigidity range, and plotted all results together. 
Data points at p < 10^'^ have been obtained from a sim- 
ulation with the full turbulence spectrum covering five 
orders of magnitude, i.e. fcmax/fcmin = 10^; the simulation 
for p > 10~^ has a turbulence spectrum that is cut-off at 
high wavenumbers but with the same ^niin '^s the previ- 
ous simulation, and fcmax/fcmin = 10'^. All rigidities have 
resonant modes with which to interact and computational 
time is saved by removing the high wavenumber modes 
(small spatial scales as compared to the Larmor radius) 
that have a negligible impact on the particle trajectory. 
This produces a slight offset in the diffusion coefficient, 
which are thus measured to ~ 10 %; however, the trend of 
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a constant D\\ / ri^c is clearly seen in both regimes p < 10^^ 
and p > 10"^. 

Note also that, in agreement with the above calcula- 
tion, the scaling of is insensitive to the slope of S^d 
in the perpendicular direction k±. On contrary, the per- 
pendicular coefRcient is sensitive to the perpendicular cas- 
cade, D± cx r^/^ . 

4.1.2. Modified Goldreich-Sridhar spectrum 

The second turbulence model we shall consider is the 
Goldreich-Sridhar model of strong MHD turbulence 
llGoldreich fc SridhadllQQ^I) modified by the streaming in- 
stability (see paper I). The three-dimensional turbulence 
spectrum then reads: 



S3dik\\,ki_) ^ So M|| ^ u^*^^ / I -1^ 



- 1 



(6) 



where f{x) is a smooth function that peaks around x = 
and vanishes at infinity. If we approximate / with a 
Dirac distribution and insert Eq.® into Eq.Q, we ob- 
tain g cx u}^^ attenuated by J^(a„) with argument a„ = 
{p\\/\n ± p±. Since the resonance condition corre- 

sponds to U|| ~ the Golreich-Shridar selection for 

the transverse modes leads to u± ^ u^^^'^ ^ '^^^ 3> 1, 
hence for short transverse wavelengths. Furthermore, for 
P\\ P± 1, the argument a„ ^ 1 and the large argu- 
ment limit of Bessel functions gives: 



9 ■ 



5/2 , 

?7P|| /PJ 



(7) 



We obtain identical results if / is approximated by a step 
function or an exponential form. 

In this anisotropic turbulence, whose spectrum is con- 
strained by the relation between k\\ and fcj^, the parallel 
diffusion coefficient D\\ shows a different scaling than in 
unconstrained turbulence. On a formal level, this differ- 
ence is related to the fact that the relation between k\\ 
and k±^ no longer permits to treat perpendicular and par- 
allel wavenumbers separately in the integrals. The QLT 
result that we obtained in Eq. Q indicates that the per- 
pendicular spectrum at short wavelengths influences the 
pitch angle scattering frequency. 

In the numerical calculation whose results are pre- 
sented in Fig. 121 we have approximated the function f{x) 
with a step function that vanishes outside of the interval 
[—1,-1-1] and used, as before, a turbulence level = 0.95. 
We find that the parallel diffusion coefficient no longer 
respects the Bohm scaling, but rather follows the perpen- 
dicular coefficient and scales as in Kolmogorov diffusion, 
D|| oc p^/^. This behavior, which does not agree with the 
above quasi- linear estimate Eq. (|7|) , may be understood in 
the following way. 

The quasi-linear calculation indicates that diffusion in 
this contrainted turbulence cannot be induced by reso- 
nant interactions in the longitudinal direction k\]^ri^ ~ 
1. Indeed, if this resonance condition is satisfied, the 
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Fig. 2. Parallel (solid dots) and perpendicular (diamonds) 
diffusion coefficients in units of tlc as a function of rigidity 

27rrL/imax- The turbulence level is 77 = 0.95 and the tur- 

5/3 1 . 

bulcnt spectrum fc^^sd oc k^ with the constraint 

that (fc^Lmax)^^'^ > fc||imax- The dotted lines have slope 

—2/3 and 1, and indicate that the diffusion coefficients 

1/3 

at low rigidities scale as as in isotropic Kolmogorov 
turbulence. 



Goldreich-Shridar prescription k\\l\\ cx: k^^^l^f^ indicates 
that fc^TL ^ 1 and this large argument of the Bessel func- 
tion in Eq. I^J then kills the contribution of these interac- 
tions in g. Since the parallel and perpendicular diffusion 
coefficients are found to satisfy D cx r\['^ , as was observed 
in unconstrained turbulence for Z?^, the numerical sim- 
ulation rather indicates that diffusion is now controled 
by perpendicular wavenumbers, i.e. k\r\^ ~ 1. If indeed 
those interactions dominate and the correlation time re- 
mains of the order of the Larmor time, the discussion in 
Appendix A shows that the observed behavior can be re- 
produced. In this description, pitch angle scattering in the 
parallel direction takes place when the particle jumps from 
one magnetic field line to another and the orientation of 
these field lines differ from one another. Diffusion along 
the parallel direction is then controled by perpendicular 
diffusion, as is apparent in the results of the numerical 
simulations shown in Fig. [51 

Following the discussion of Appendix A, we finally ob- 
tain: 



VP 



2/3 



(8) 



It is interesting to note that the exponent does not de- 
pend on the turbulence spectral indices a and /3 in a first 
approximation (see Appendix A). 

For > 1, the arguments of the Bessel function are 
small and we find again the same result as in the isotropic 
case, namely a parallel diffusion coefRcient increasing like 
p|. Interestingly, the above diffusion coefficients should 
be expected to apply to propagation in the interstellar 
medium as well if the streaming instability of CR along 
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the mean magnetic field is the dominant source of turbu- 
lence generation. A detailed investigation of this issue is 
postponed to a future work. 



4.2. Downstream diffusion coefficients 

Downstream of the shock the parallel and perpendicu- 
lar coherence lengths are different: they are related to 
the coherence length upstream Lmax by £_l = ^max and 
^11 = Lmax/?'- The turbulence transforms in a non-trivial 
way, since the magnetic field modes components that are 
perpendicular to the shock normal are amplified by r, and 
these modes may be associated with wavenumbers that 
are neither completely parallel nor completely perpendic- 
ular. The overall turbulence magnetic field is nevertheless 
amplified by i? = ^(1 -t- 2r'^)/2> after statistical averaging, 
while the coherent magnetic field, which was assumed par- 
allel, is left unmodified. Hence the turbulence amplitude rj 
changes from upstream to downstream by the above am- 
plification factor Rii rj <^ 1, or remains the same if 77 ~ 1. 
The effective value of the compression ratio explored by 
the particles depends on their rigidity in the case of a 
strongly modified shock. Only the CRs of highest energy 
should explore the total compression ratio, rtot while low- 
energy particles explore the much lower compression ra- 
tio of the subshock, rgub, and the electrons of highest en- 
ergy explore a region across the shock with an interme- 
diate effective compression ratio jvi- The fluid quantities 
(density, velocity, thermal pressure, tangential magnetic 
field) are all subject to a jump dependent on rgub at the 
shock precursor. For instance the magnetic field compres- 
sion R ratio should be expressed in terms of Tgub instead 
of r. Typical values are rsnh 9^ 2 — 3 and rtot — 7 — 10 
l|Berezhko et aLlllQflfl iBvkov 1120041) . However, in section 
El we will show that a consistent modelling of young SNr 
with a strong magnetic field amplification by the stream- 
ing of CR should include both the CR backreaction on the 
shock structure and the CR distribution steepening pro- 
duced while generating the turbulence. The last effect will 
tend to diminish the CR pressure and to reduce the non- 
linear shock modification. One may expect in that case 
Tsub not to be strongly different from the the test particle 
case. 

Mainly for the previous reason and also as they would 
only complicate the formal treatment with no significant 
change in the results we shall not go into such details 
in this paper. (Accounting for a different compression ra- 
tio Tsub and Tel would lead to a variation of the order of 
30 % for the deduced v alue of i?p,max as already stressed 
in lParizot et al ](l2006*).) In any case, we found that a sit- 
uation with explicitly different values of the compression 
ratio for nuclei and electrons in the case of a strongly mod- 
ified shock results in an intermediate solution between the 
test-particle case (no precursor nor sub-shock, and r = 4) 
and the pure r = 10 case e xpected in strongly modi- 
fied shocks fsee lBvkovl l|2004l) 'l. We shall thus confine the 
study below to such idealised cases, using only one ("uni- 
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Fig. 3. Same as Fig.^ i.e. a turbulence spectrum k^S^d oc 
kj^^^^k'^^ and 77 — 0.95 but accounting for the compression 
of the turbulence through the shock. 
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Fig. 4. Same as Fig. |21 i.e. assuming that (fc^imax)^^"^ > 
^ll^max, and accounting for the compression of the turbu- 
lence through the shock. 



versal" ) compression ratio r for all particles and the fluid 
quantities. 

The transport properties differ in co mpressed and uncom- 
pressed turbulence, see for instance iLemoine fc Revenu I 
(200.5.) for the case of ultra-relativistic shocks. In the fol- 
lowing Figs. l3l4l we show the diffusion coefficients of parti- 
cles propagating in compressed turbulence (with r = 4) for 
the same cases of turbulence spectra as shown in Figs. lll2l 
The comparison of these figures reveals that the main 
effect of compressed turbulence comes through the reduc- 
tion of the parallel coherence length: in compressed turbu- 
lence, the scaling of Z?|| with respect to rigidity 27rrL/Lfnax 
remains the same as in uncompressed turbulence, but the 
curve is shifted to lower values of the rigidity. Hence if one 
plotted Z^ii versus r x 27rrL/imax, most of the discrepancy 
would disappear. One also observes that the parallel, and 
to some a smaller extent the perpendicular diffusion co- 
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efficients are decreased to a lower value. The difference 
in turbulence level between downstream and upstream is 
probably too small to account for tliis difference, since 
7] = 0.95 upstream and r] = 0.98 downstream. In tfie case 
of a modified Goldreich-Sridhar spectrum, a reduction fac- 
tor of ~ 42/3 of the maximum turbulence length can ex- 
plain this difference. 

Another consequence of magnetic field amplification 
and compression at the shock front is that the ratio of the 
downstream Alfven velocity VacJ to the downstream fluid 
velocity Vd is unlikely to be much smaller than unity as in 
the upstream medium far from the shock front. In detail, 
the downstream Alfvenic Mach number A^Ad = Ki/VAd is 
expressed in terms of the upstream Alfvenic Mach number 
Mau as: 



(9) 



due to compression of the magnetic field and density by 
~ r. With r ~ 4 and A^au ~ 10 [see Eq. (^Hl below], one 
finds VacI ~ Ki. In this case, second-order Fermi accelera- 
tion should be taken into account during the Fermi cycle. 
The stochastic momentum diffusion coefficient is easily 
expressed in the relativistic limit as 



2 / Va 



9 ^ Dn 



(10) 



and p — e/c is the particle momentum. The Fermi sec- 
ond order acceleration timescale, tpii = [Dpp/p^)^^ oc D\\ 
is then controled by the downstream stochastic diffusion 
coefficient. This contribution of stochastic reacceleration 
downstream is discussed in Sectional 

5. Revised Fermi cycle in shocks with magnetic 
field amplification 

Each relativistic particle streaming upstream of the shock 
front generates shear Alfven waves. Accounting for the 
possibility of strong turbulence, the energy loss in gener- 
ating the waves is 



P{e) 



1 

'3^ 



V 



d_ 

dx 



log/ 



(11) 



(see Appendix B), V the velocity of the scattering centers 
upstream is close to the local Alfven velocity Vau- 

This general expression can be specified upstream 
where the distribution function is governed by the follow- 
ing transport equation, assuming a stationary state: 



d d d 
ox ox ox 



(12) 



The power loss, in most of the upstream region where 
f ~> fvL, is then given by 



(13) 



The average duration of a half Fermi cycle upstream is 
4£)u/ (Vuw). A particle that undergoes such wandering and 



thus contributes to turbulence generation spends an aver- 
age amount of energy — The energy gain per Fermi 
cycle G = 1 -I- ^ is thus diminished for two reasons, first 
by the fact that the particles experience a reduced veloc- 
ity difference between scattering centers V^i — V^ — V and 
moreover loose energy by exciting turbulence; then 



4K_ 
3 



Vd 



Vv 



(14) 



Since the escape probability per cycle is unaffected, 
-Pcsc = 4Vd/f, the amplification of upstream turbulence 
and the concomittant reduction of the energy gain steepen 
the accelera ted spectrum. The spectral index is given by 
(BdOElzi): 



^ ^ ^ log l/Prct _ ^ ^ Pcsc_ _ ^ 



logG 



s ~ 1 



[r - 1) Vsh 



(r-1) 



, (15) 



(16) 



or again 



2,8 Va^ 
3 ^sh 



(17) 



for a compression ratio r = 4. From Eq. (|15|l . we can see 
that the effect of wave generation upstream on the CR 
spectrum is important once the amplitude of the fluctu- 
ations 5B becomes much larger than the ISM magnetic 
field. Let us quantify this remark. The index can be ex- 
pressed in terms of the inverse of the turbulent (upstream) 
Alfven Mach number 



Ma 



14.3 



B 



lOO^iG 



1 cm 



1/2 l^h 



10- 



.(18) 



Instead of an index of 2 for strong shocks, we get 
for a SNr shock Vsh — 5000 km/s, with an amplified 
upstr e am magnetic f i eld of ~ 100/^G ijBerezhko fc Vol^ 
l2004t IParizot et all l2006|) in a hot ISM medium 
riHiM — 4 10~3 cm~'^, s = 2 + 1.5 and in a warm ISM 
medium nwiM — 10^^ cni^"^, s = 2 4- 0.35 respectively. 
The index appears to be strongly modified especially in 
the hot phase. 

The compression ratio can also be changed with a cor- 
rection of order 1/A1au by the modification of the shock 
by both CR-pressure gradient and wave-pressure gradi- 
ent. Note that the thermal plasma pressure also increases 
in the precursor because the magnetic energy cascade to- 
wards dissipation scale where cyclotron absorption takes 
place. These effects tend to diminish the compression ra- 
tio. Using Eq. (|18|l for r = 10 the spectral index s ~ 1.65 
in the warm ISM is found to be slightly harder than the 
solution obtained in modified shock simulations (s ~ 1.5, 
lEUi son et al. I ('2004')). All these effects arc still to be in- 
tegrated in numerical simulations of CR modified shocks 
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including kinetics an d the effect of a strong Alfven pro- 
duction upstream [see lJone s (1994) for a discussion in the 
three fluid approximation]. 

As already mentioned, the generation of a highly am- 
plified and compressed magnetic field enhances the Alfven 
velocity compared to far upstream in such a way that 
the Fermi second ordrer process may become important. 
However, if only forward waves are generated, the Fermi 
second order process does not work. Assuming that the en- 
ergy is distributed evenly between forward and backward 
waves, one can estimate the mean energy gain through 
second order acceleration using Eq. ll()|l : 



5FII 



dp 



(19) 



The overall particle energy gain/loss during one cycle 
transforms to: 



Gtot — 1 + T 



4K-Vd 4 [Vau , VauV 



.9FII 



(20) 



If we note gpi = 4/3 (Vu — Vd)/v the mean energy gain 
corresponding to the first Fermi process, then the use of 
Eq. (^nil gives: 



gFII 
9F1 



3 — K 
r - 1 



(21) 



where k is defined by the energy (rigidity) dependence 
of the diffusion coefficient: Dy cx -DBohmP'' oc p*^"*"^, with 
-DBohm = ''lc. If the anisotropy of the turbulence is not 
constrained by a relation of the form fcii dis- 
cussed in Section [4. 1.11 then k = 1 — /?, with j3 the index 
of the turbulence power spectrum as a function of fc|| . In 
particular, (3=1 implies k = 0, hence a Bohm diffusion 
coefficient. 

Usually, magnetic field strength up to 500 /iG are de- 
rived in young SNr producing non-thermal X-ray radia- 
tion. The previous equation implies that the second order 
process can become comparable or even dominate the reg- 
ular Fermi process for ISM densities lower than a few tens 
(corresponding to the hot ISM phase) for such high mag- 
netic field strengths. The stochastic process has then to 
be incorporated in a self-consistent way in the Fermi cycle 
in such cases. However, in order to do so one has to take 
into account the likely anisotropy of the turbulence be- 
tween backward and forward waves and evaluate its effect 
on the efhciency of the second order Fermi process. This 
investigation remains beyond the scope of this paper and 
is left for future work. 



6. Discussion: Cosmic-ray acceleration in SIMr 

In this section, we aim at deriving the maximum CR ener- 
gies produced by the diffusive shock acceleration processes 
in an anisotropic turbulence. We limit our calculations to 
the young SNr showing thin X-ray filaments as observed 
by the satellites Chandra and XMM-Newton (see iBalleti 



iooi)). We postpone to a future work the derivation of 
the overall source CR spectrum produced by SNr during 
their whole evo lution phases as discussed recen tly in a sim- 
ilar context bv lPtuskin fc Zirakshvilh I l|200fj) . 
In section 4, we have derived the general form of the diffu- 
sion coefficient up- and down-stream. Expressed in terms 
of the Bohm diffusion coefficient i^Bohm = ^hC, we get: 



D(p) = fcDBohm — ,forp<pp 



= fcD 



Bohm 



, for 



P> Pp 



(22) 



From the simulations we infer approximatively ppu — 
0.2 upstream and — 0.2/r downstream, r being the 
shock compression ratio. The normalisation factor k is 
also different up- and downstream. In the case of a Bohm 
type turbulence upstream, compressed by a factor r down- 
stream (assumed to be 4 hereafter) we have ku ~ 2 and 
kd — 4 respectively (see figures 1 and 3). In case of a 
Kolmogorov type turbulence upstream, compressed by a 
factor r downstream, k^ ~ 6 and kd ~ 6/4^/'^ ~ 2.4 (see 
figures 2 and 4). 

In order to discuss the particle acceleration process in 
young SNr quantitatively, all the lengths are compared 
to the size of the SNr and all the timescales to the age of 
the SNr. The maximum scale of the turbulence up-stream 
and the particle gyroradius can be expressed in terms of 
the SNr shock radius units, namely Lmax = Lmax x R-sh 
fL = fL(E,B) X Rsh with, fL = fo Epcv ^ioq, Epev is 
the CR energy in PeV and Bioo is the magnetic field in- 
tensity in 100 /iCauss units. Hereafter, we will restrict 
our analysis to young SNr in a free expansion phase for 
which the forward shock velocity is the largest. This leads 
to a direct relationship between the SNr radius and age 

Rsh — Vsh X tsNR,: 



R 



2^^ 



sh 



tSNR 



5000 km/s 300yr ' ^'^^^ 

hence f o — 2 10~^ for these values. We now derive the 
explicit values of the maximum cosmic-ray energies. 

Note that downstream, the turbulence m ay also re- 
lax t owards ISM values over a scale ii-(E) l|Pohl et al I 
l2005l) producing a diffusion coefficient Dd — D^^ir 
oo) X exp(x(a-|-l)/2£r), with Dd{ir — ^ oo) the diffusion co- 
efficient in a uniform medium as before. Relaxation scales 
^ Rsh have for average effect to increase the particle 
residence downstream and then the acceleration timescale. 
It is expected that acceleration through a relaxed turbu- 
lence is less efficient in producing high energy CR. A de- 
tailed investigation of this process deserves to be pushed 
further. 

Maximum CR energy in young SNr: 

If the relaxation scale is ~ Rsh then the X-ray filaments 
observed in some young SNr are limite d by the radia- 
tive (synchrotron) losses. Several authors jBerezhko et alJ 



l2003HBerezhko fc V61kl2004l:IParizot et al 12006^ have de 
rived in a self-consistent manner the downstream mag- 
netic field by comparing the observed filament size Aiiri„ 
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with combined advection and diffusive lengths explored 
by the relativistic electrons during their synchrotron loss 
timescale. The magnetic field strengths obtained are of 
the order of 300 — 400 /zG in the youngest SNr like Kepler 
and Tycho and 100 /LtG in older SNr like SN1006. These 
estimates however have assumed an isotropic turbulence 
up- and downstream. 

As discussed in the introduction, the anisotropy 
considered here has two main effects. Firstly, if the 
Goldreich-Sridhar scaling does apply then diffusion has a 
Kolmogorov scaling even if the streaming instability tends 
to produce a spectrum oc ^ij"^. Secondly, this turbulence 
spectrum is transmitted downstream and the scales par- 
allel to the shock normal (in a quasi-parallel shock config- 
uration) are compressed by a factor close to r. We shall 
account for both these effects in the following estimates. 

We first derive the magnetic field downstream us- 
ing th e size of the X-ray fi l ament s [see iBerezhko fc Volkl 
l)2004|) and IParizot et all l)2006|) ]. The electron parti- 
cle distribution f{x) verifies Vd{df/dx) ~ Dd'^f/dx^ — 



fh 



In this equation, catastrophic losses for TeV elec- 



trons have been assumed. The solution is of the form 
f{x) oc exp(— x/Ai?i.i„i). The magnetic field downstream 
then follows from: 



D 



Ai?riT 



1 



'syn 



(24) 



The synchrotron loss time r^yn ^ (95yr) B^^^^^ ^2l-\^e\^ 
where i?obs-kcV is the photon energy in keV units at which 
the X-ray filaments are observed and Bdioo is the down- 
stream magnetic field in units of 100 /iG. Once the de- 
pendence of the diffusion coefficient on particle energy is 
known, Eq. H24|l leads to a one to one relationship between 
the downstream magnetic field and Ai^rim = Ai?rim ^sh- 
The ratio of the upstream to downstream magnetic field is 
approximately 1.2/r (Parizot et al 2006). The maximum 
cosmic-ray energy -Epmax is then calculated by balancing 
the diffusive length upstream D-alVs\x with the shock ra- 
dius [see Hillas (20Q5|) for a discussion]. In order to fix 
the maximum scale of the turbulence, let us first assume 



Bohm regime: In the case where diffusion proceeds accord- 
ing to a Bohm regime upstream, k = 0. The magnetic 
field downstream is then derived as ~ 400 /iG for a 
typical filament size Ai?i-i,„ = 10~^ and r = 4. The results 
are similar for r=10 in case of a strongly modified shock. 
Such high magnetic fields downstream require an efficient 
amplification in the shock precursor (see paper I). Using 
the above argument on upstream CR escape, we find max- 
imum cosmic ray energies i?pmax — 1.5 PeV for r = 4 or 
0.6 PeV for r 10. 

Kolmogorov regime: if the upstream turbulence is of 
the Goldreih-Shridar type, the diffusion is similar to 
Kolmogorov, k, = —2/3 and Ea. H24|l leads to a mag- 
netic field of the order of ~ 550 /zG for filament size 



Ai?Hm = 10"^ and r = 4, Bd ~ 400 /iG for r = 10. 
Using the above argument on upstream CR escape, we 
find maximum cosmic ray energies -Epmax — 7 PeV (r — 4) 
or ~ 2PeV (r = 10). 



These estimates suggest that it is possible to accelerate 
CR up to the CR knee with magnetic field amplification 
in the shock precursor for both Bohm and Kolmogorov 
regimes. Th is confirms earlier r esults found for isotropic 
turbulence l)Parizot et all 1200^ even if the typical max- 
imum energies for CR tend to be slightly higher in case 
of anisotropic compressed turbulence. Another difference 
with the isotropic turbulence case, is that the Kolmogorov 
regime can not be ruled-out here; the rejection condition 
D < i^Bohm is not true anymore in anisotropic turbu- 
lence. However, it seems still challenging at this SN evo- 
lution stage to accelerate the particles up to t he CR an- 
kle [see Ptuskin & Z irakshvilli (2005), Bvkovl l|200lh and 
IParizot et al r|2004l) for alternative scenario]. 

The above maximum CR energies are probably slightly 
overestimated since the diffusion coefficient oc at 
high rigidities and the hypothesis Lmax — ?'L(£'pmax) is 
optimistic. In order to estimate this uncertainty, we use a 
supplementary relation to fix imax- For instance, the con- 
dition racc(i?cmax) = Tgync (£^emax) at thc maximum elec- 
tron energy provides a one to one relationship between 
^max and the magnetic field. For the youngest SNr, the 
typical synchrotron cut-off frequency is of the order of 
1 keV. Downstream magnetic fields in both regimes are 
still of the order of 400 — 500 /iG. In the Kolmogorov 
regime, Lmax — 10~^ and the maximum CR energies 
cannot lie well beyond 0.3 PeV as a result of the scaling 
Z) oc at high rigidities. 

To summarize, the above estimates give as an order of 
magnitude for the maximum energy of cosmic rays: E ~ 
Z X (0.3 — 3) PeV in young SNr showing X-ray filaments if 
the magnetic field is amplified by the streaming instability. 

This paper only considered the free expansion stage 
of SNr evolution and did not investigate the overall 
spectrum, expected to be st eep beyond the CR knee 
l)Ptuskin fc Zirakshvilli1l2005l) . The maximum CRs ener- 
gies calculated here are estimates expected once a strong 
magnetic field amplification is at work at the shock precu- 
sor as deduced from X-ray observations from young SNr. 
Further work should take into account the different SNr 
evolution phases either, as well as the different type of SN 
(core-collapsed or type la). A consistent spectrum deriva- 
tion should account for both upstream wave generation 
and shock smoothing by the cosmic-ray pressure as well as 
the second order Fermi acceleration process downstream. 
Consequently, a detailed reconstruction of the CR spec- 
trum produced by an ensemble of SNr during their evolu- 
tion is beyond the scope of the paper and will be adressed 
in future work. 
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7. Conclusion 

In this paper II, we have shown that the diffusion coefh- 
cients of cosmic rays can be significantly modified in the 
precursor of a shock, especially because of the formation 
of an anisotropic turbulence spectrum. According to paper 

1, the streaming instability shapes the power spectrum of 
the turbulence according to the scahng S^d cx k^^kj^"^^, 
while the transverse index a depends on the modeling 
of the Alfven cascade. In the absence of any correla- 
tion between parallel and transverse wave numbers, the 
spectrum factorizes and the diffusion coefficients are in- 
sensitive to the transverse distribution. The parallel dif- 
fusion coefficient respects a Bohm scaling, D cx tlc. 
However, if a correlation a la Goldreich-Shridar arises, i.e. 
fc||£|| ~ (fc±^j.)^^'^, the results differ from the quasi-linear 
prediction. Numerical simulations indicate that the paral- 
lel diffusion coefficient scales with rigidity as in isotropic 
Kolmogorov turbulence. This behavior is understood if 
pitch angle scattering is now controled by transverse trans- 
port from one field line to another. In any case these re- 
sults allow to explain the recent X-ray observations young 
Supernova remnants. The diffusion coefficient for the elec- 
trons, in particular, agrees well with the expected size of 
the filaments. 

The energy loss per particle that accompanies the gen- 
eration of turbulence is calculated. The cosmic-ray en- 
ergy spectrum is consequently steepened, which supports 
the suggestion that the index at the source is closer to 
2.3 (index obtained in the warm ISM phase) rather than 

2. A further steepening due to the diffusive propagation 
in the interstellar medium would likely leads to an in- 
dex close to 2.7 as long as the diffusion develops with 
a Kolmogorov spectrum. Finally the amplification of the 
magnetic field allows to push the maximal energy cut-off 
into the "knee" region (PeV energies). Consistent future 
investigations with the cosmic ray data should require the 
inclusion of the magnetic field amplification by the stream- 
ing instability, the shock modification by the cosmic ray 
pressure and possibly the effect of stochastic Fermi accel- 
eration downstream the shock. 
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Appendix A: Pitch angle diffusion 

The variations of the pitch angle are governed by a simple 
stochastic equation that stems from the projection of the 
Lorentz equation along the mean field by taking into ac- 
count the energy conservation (and thus the conservation 
of p and v). This gives: 



where lol = ZeB /rwyc^ (j){t) is the gyro-phase, i.e. (j){t) — 
lol + 0{h) and b = SB/B is the irregularity of the field 
experienced by the particle along its trajectory. For a sta- 
tionary process, the pitch angle frequency is thus 



At 



2 1 (/(r)/(0))dr , 

'0 



which can be rewritten as follows: 

/>OC 

i^s=(^l / (b(r) • b(0) cos A0(t)) dr 



(A.2) 



(A.3) 



If the level of turbulence is sufficiently weak, one approx- 
imates A0(r) ~ w^T and b(T) is expressed in term of the 
unperturbed trajectory. This leads to the expression: 



d^k 
(2^ 



S3dik)Rik,p) 



(A.4) 



where 5*3(1 (fc) is the 3D correlation spectrum of the field 
irregularities normalized such that its integral equals to 
the irregularity degree rj: 



{SB^ 



(Pk 
(2^ 



(A.5) 



and the resonance function R{k,p) describes the resonant 
interaction between the particles and the modes: 



r 

R{k,p) = ujl I 
Jo 



^ik.Ax{T)—iiij{k)7 



cos(wiT) dr , (A.6) 



where Ax{t) is the variation of the unperturbed trajec- 
tory during a time lapse r and w(fc) the mode pulsation. 
This exhibits the Landau-synchrotron resonances of the 
form: 

R{k,p) oc S[uj{k) - ki\vi\ ±nuJL] ■ (A. 7) 

For an isotropic power law spectrum S{k) cx rjk^^^ , 



(A. 



where p = 27rrL/imax, ^max being the coherence length 
of the field. The parameter p is the reduced rigidity and 
this law holds for: 



jO ^ 1 



(A.9) 



Jit) = OJL [cOS(P{t)b2{t) ~ Sm(j){t)bi{t)] , 



(A.l) 



Finally the scattering time Tg ~ t's"^- 

In the case of strong turbulence, the gyro-resonances 
broaden; it turns out that the scaling law in terms of the 
rigidit y p and of the irre gularity level rj can be extrapo- 
lated 1 Casse et al.l l|2nn2l) 'l when the spectrum is isotropic. 
In this paper, we show that the results can be different 
when the spectrum is anisotropic. When there is no corre- 
lation between the parallel and the perpendicular parts of 
the spectrum (i.e. S'3d(fc) cx k'^^kj^ with q = a + 1 > 2) 
the result derived from the quasi-linear theory can be ex- 
trapolated to strong turbulence, giving a diffusion law sim- 
ilar to the isotropic case. In particular a Bohm scaling is 
found for /3 = 1 (and only in that case). 
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The result is quite different when there is a correla- 
tion between parallel and transverse wave numbers as pre- 
scribed by Goldreich and Shridar, i.e. k±£± ^ 



fnserting ea. l|B.3|l and ca. (|B.2p into ca. (|B.ip . one obtains 



Now the resonance fc||rL ~ 1 becomes inefficient because / P(e)/(p)47rp^dp = / (1 — fi^)d^ 

the requirement to have simultaneouslv much larger k\ J ^ J-i 



the requirement to have simultaneously much larger fcj 
makes the particle insensitive to these modes (i.e. the an- 
gular scattering frequency vanishes). The numerical sim- 
ulation actually indicates that the particles interact at 
fcj^TL ~ 1, as shown in Fig. |21 This behavior is not de- 
scribed by the quasi-linear theory that would give a non 
vanishing result only if one could obtain simultaneously 

/CIITL ~ f. 

However, if one assumes that the particle interacts 
with a spectrum band such that k±ri^ ^ 1 and the corre- 
lation time is a few Larmor times, as is usually the case, 
Eq. (jSU gives a relation of the form: 



dk 



(2^) 



S3d{k±, fci 



(A.fO) 



fcj_l'L = l 



which leads to a result in agreement with the numerical 
simulation, if one recalls t he rel ation 3a + 2/3 = 7 [see 
paper I and lCaltier et al. I ||200,'tI) ]: 



2/3 



(A.ll) 



This effect can be associated with a perpendicular diffu- 
sion process associated to encounters with oblique mag- 
netic field lines encounterings. The particle angular diffu- 
sion is then controled by the perpendicular transport. 

Appendix B: Power loss rate per particle 

Let 7r(fc) be the growth rate of the resonant instability 
(see paper I) and P(e) the power lost by each particle of 
energy e to trigger the turbulence spectrum Sad (A;). Then, 
these physical quantities are linked the following integrals, 
expressing that the power generation of waves is equal to 
the power loss suffered by the particles: 



P(e)/(p)4Vdp = -2 / 7,(fe)53d(fc) 



d^k 
(2^ 



.(B.f) 



The growth rate is proportional to the gradient of the 
distribution function and depends on the pitch angle fre- 
quency Vs(p,^). As seen in section 0] the spectrum fac- 
torizes (except in the Golreich-Shridar case) so that the 
scattering frequency is independent of the transverse spec- 
trum and is given by 



z/s(p,Ai) = ^i.{v)k\\^\\{k\\)\ 



(B.2) 



where fci.(p, /i) = (|/i|rL) ^ ■ Then the growth rate is given 
by the following integral: 



■kV 



7r(fc||) = ^2 P (1 - M 

^,2 Qj- 

x_e5(fc||rLM-l)^. 



(B.3) 



ox 



(B.4) 



This result is in agreement with ijMcKenzie fc Volkl 
if 982)) who obtained a power density for Alfven wave gen- 
eration equal to VdPcr/dx. One finally obtains the power 
loss per particle: 



(B.5) 



However a model of CR acceleration at SNr shock to 
be consistent with the CR data should include at once 
the amplification of turbulence upstream and the back 
reaction of cosmic ray on the shock structure. 
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